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Abstract. We formulate gaussian and circular random-matrix models representing 
a coupled system consisting of an absorbing and an amplifying resonator, which are 
mutually related by a generalized time-reversal symmetry. Motivated by optical 
realizations of such systems we consider a VT or a VTT' time-reversal symmetry, 
' (-h '. which impose different constraints on magneto-optical effects, and then focus on five 

| common settings. For each of these, we determine the eigenvalue distribution in the 

complex plane in the short- wavelength limit, which reveals that the fraction of real 
eigenvalues among all eigenvalues in the spectrum vanishes if all classical scales are 
kept fixed. Numerically, we find that the transition from real to complex eigenvalues in 
the various ensembles display a different dependence on the coupling strength between 
the two resonators. These differences can be linked to the level spacing statistics in 
the hermitian limit of the considered models. 
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1. Introduction 

The investigation of nonhermitian "PT-symmetric Hamiltonians is motivated by the fact 
that they possess eigenvalues which are either real or occur in complex-conjugate pairs 
p] . Considerable attention has been paid to the delineation of systems with a completely 
real spectrum, with many works focussing on exactly solvable one- dimensional situations 
(for reviews see [2j [3]). With the recent advent of optical implementations jH [5] it has 
been realized that the appearance of complex eigenvalues drives a number of interesting 
switching effects [U El [7J El [9J [TUl HH [12], including the possible onset of lasing 
[T3l dH [151 HE] , which moves the most unstable states (with energies or frequencies that 
have a large positive imaginary part) into the centre of attention. At the same time, 
these implementations motivate the study of multi-dimensional systems in which many 
modes become mixed by multiple scattering. Here, we investigate the formation and 
distribution of the complex spectrum in such situations on the basis of a statistical 
approach rooted in random-matrix theory [TT1 IT8] , which samples systems that respect 
a certain set of symmetries and share a number of well-defined characteristic energy and 
time scales, but differ in the microscopic details of the dynamics. 

We extend earlier exploratory works of this approach [191 [20J to consider random- 
matrix ensembles which differ by the assumed absence or presence of elastic or 
dissipative magneto-optical effects. This leads to a choice between two generalized 
time-reversal symmetries, termed VT and VTT' symmetry and physically motivated 
in [2T]. These ensembles apply to a coupled-resonator geometry (with an absorbing 
resonator possessing M internal modes coupled to a matching amplifying resonator 
via an interface of N channels with transparency T, and amplification or absorption 
rate set to a common value /i). The optical setting motivates to consider 5 particular 
scenarios (00, UO, UO', OA and OA'). These can be studied either based on an 
effective Hamiltonian or in terms on an effective time-evolution operator (a quantum 
map), as is described in section [2j 

In section |3] we determine for each scenario the distribution of eigenvalues in the 
complex plane in the short- wavelength limit M — > oo at fixed a = M/N, T and fi. We 
find that in this limit, the fraction of real eigenvalues among all eigenvalues vanishes 
at any finite fixed amplification and absorption rate, with the details of the eigenvalue 
distribution in the complex plane depending on the symmetry class. This supports the 
conclusion of earlier work on some of these ensembles [19] that the transition to the 
complex spectrum occurs at a characteristic absorption rate /ipx which is classically 
small when compared to the inverse dwell time E<i = 1/tdweii = NTA/27T in each part 
of the resonator, but large when compared to the mean level spacing A. 

In order to investigate the details of this transition we then present results of 
extensive numerical investigations (section H] and [5]). These reveal that the various 
ensembles display a different dependence of the transition on the coupling strength 
T, as well as on M and N. This division is associated with specific mechanisms of 
eigenvalue coalescence, which we relate to differences in the level spacing statistics in 
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Figure 1. Illustration of (a) the effective Hamiltonian ([1) and (b) the quantum map 
([6]) used to model an absorbing resonator (L) which is symmetrically coupled (via an 
interface characterized by T or C) to an amplifying resonator (R). Various symmetry 
classes arise depending on the constraints imposed on the internal Hamiltonians 
and The depicted situation applies to uniform amplification or absorption with 
rate [i; further symmetry classes arise when /i is replaced by matrices Xl and X-&. 



the hermitian limit by extending the perturbative considerations of [19J. 
Section [6] contains our conclusions. 

Throughout this work we denote eigenvalues as E, but set H = 1; all considerations 
thus directly apply to the eigenfrequencies in optical analogues of non-hermit ian 
quantum systems. 



2. Random-matrix ensembles 

Following [19j [201 EI] , we consider a coupled- resonator geometry where one part of the 
system (L) is absorbing and the other part (R) is amplifying, with the absorption and 
amplification rate set to a matching value [i. In each part random multiple scattering at 
a rate 1/r results in a mixing of M = 1/Ar internal modes, where A is the mean level 
spacing, and the two parts are coupled together at an interface which supports N open 
channels of transparency T n , n = 1, 2, . . . , N. In order to capture the consequences of 
multiple scattering we utilize effective Hamiltonians and quantum maps, which model 
these systems as illustrated in figure [H 



2.1. Effective Hamiltonians and symmetry classes 

The general structure of an effective Hamiltonian for this situation can been derived in 
a systematic scattering approach [191 EI] • This yields 

where the M x M-dimensional hermitian matrices and (Xl and Xr) represent 
the internal hermitian (anti-hermitian) dynamics in each part of the system, while the 
coupling matrix is of the specific form 

AM y/T^ 

N entries M — N entries 

The specific form of ([T]) displays a structure which goes beyond the mere symmetry 
requirements usually applied in mathematical classifications of nonhermitian matrices 
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(for comprehensive overviews see [2"2"| I2"3"t 124]). In particular, the matrix T needs to 
be positive definite and bounded in order to model physical coupling between two 
resonators. This structure resembles analogous models in mesoscopic superconductivity, 
where the two subspaces represent particles and holes, and the coupling is provided by 
Andreev reflection [251 [26| [27] . 

We now impose two different versions of generalized time-reversal symmetry. 
Traditional VT symmetry involves the parity operator V = o~ x ® 1m, where the Pauli 
matrix o~ x interchanges the subspaces R and L, as well as the time reversal operation 
T = K, where K is the complex conjugation in a given basis, assumed to coincide with 
the basis of ([I]). Invariance under the joint VT operation then demands 

H = VH*V =► H h = H* K = Hl, X L = X*=X£. (3) 

In a "PT-symmetric basis, the secular polynomial s(E) = det (H — E1 2 m) has real 
coefficients, which constraints each eigenvalue to be real or being partnered by its 
complex conjugate. 

In hermitian situations, the complex conjugation T is equivalent to taking the 
transpose of the matrix (thus passing from the right eigenvalue problem to the left 
eigenvalue problem). In non- hermitian situations, this transposition amounts to an 
independent operation, denoted as T' [191 EI]- For a 'PTT'-symmetric situation, we 
now obtain the constraints 

U = VV)V => H L = H R , X L = X R} (4) 

which is of interest as this yields the same spectral constraints as VT symmetry. 

For each of these two number of ensembles can now be formulated depending 

on the presence or absence of additional symmetries for H = and X = X^. In 
particular, we consider the cases that they may be further constrained to be real and 
thus symmetric (labeled O for orthogonal), complex (labeled U for unitary), or purely 
imaginary and thus antisymmetric (labeled A). In combination, we then arrive at 9 
symmetry classes with "PT-symmetry, denoted as S^Sx, Si = 0,U,A, as well as 8 
additional classes SfiSx with "PTT'-symmetry (00 and 00' coincide as in this case T' 
is an independent symmetry). 

A detailed discussion of the physical requirements corresponding to the various 
symmetries in optical settings is given in [2]]. Motivated by this context, we focus 
on 5 situations, 00, UO, UO', OA and OA'. The most important scenario is that of 
00 symmetry, with H = H T = H*, X = X T = X*, which includes ordinary optical 
systems with gain and loss modeled by a complex refractive index. The cases of UO 
and UO' symmetry (X = X T = X* but H not further constrained) model systems 
with elastic magneto-optical effects, with different symmetry constraints imposed on 
the effective magnetic field. We will also consider the OA and OA' it is known 

that absorption can be provided by magneto-optical devices [28] (in practice, however 
the design of a matching magneto-optical amplification may prove challenging). 
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2.2. Random-matrix ensembles 

The described symmetry classes are converted into statistical ensembles under 
convenient sampling of the M x M-dimensional hermitian matrices H and X. 
Specifically, depending on whether Sn = O or U we choose H from the standard 
Gaussian orthogonal or unitary ensemble (GOE or GUE) of random matrix theory 
[171 EE], respectively. The variance \Hi m \ 2 = a = 1/M of the matrix elements is set such 
that the probability distribution of eigenvalues E becomes stationary in the large-M 
limit, corresponding to a Wigner semicircle law with radius 2, 

p{E) = 7T- Vl " E 2 /A. (5) 

For Sx = O symmetry of the anti- hermitian part, we model uniform absorption 
and amplification by setting X = [iIm- For the case Sx = A, we model iX = —A via a 
real antisymmetric matrix with random Gaussian elements, and quantify the degree of 
non-hermiticity by fi 2 = M~HxAA T '. 

Throughout, we will model the coupling between the two parts of the system via 
N = aM channels of identical transparency T. Together with the chosen energy scaling 
(jSJ), which gives Ao = MjpiE = 0) = M/tt, the coupling matrix (j2j) then takes the 
form T = diag(7, . . . , 7, 0, . . . , 0) = 7 diag(l7v, Om-n), with N finite diagonal entries 

7 = VT(i + y/T=T)-\ 

We denote these ensembles as GS'h-S'xE or G^h^xE', and specifically consider the 
cases GOOE, GUOE, GUOE', GAOE, and GAOE', which correspond to the optical 
settings described in the previous subsection. 

2.3. Effective quantum maps 

An alternative approach in random-matrix theory bases the considerations on circular 
ensembles of effective time-evolution operators [TTt IT8] . For the coupled- resonator 
geometry, the general structure of these operators has been identified in [2"0j l2"Tj . They 
take the form of a quantum map 

r-( e~^E L \ r- r- ( Re^P + Q -iIm7P \ 
^=^1^ e-F R J^' V °=[ -ilmjP Re7P + gJ' (6) 

which delivers quasienergies E n via the eigenvalue problem 

Tip n = X n ip n , X n = exp(-iE n r). (7) 

The properties of the interface are now encoded in the parameter 7 = a/ y/R + iVT, the 
rank-A^ projector P = diag(l7v, Om-n), and the complementary projector Q = 1 M — P. 
The internal dynamics are described by the M x M-dimensional unitary matrices 
and Fr, which satisfy F^ = F£ for VT symmetry, and F^ = for VTT' symmetry. 
Finite \i breaks the unitarity of the quantum map. (In the specified form, (jfJJ) holds for 
Sx = O symmetry with uniform amplification and absorption, but by the replacement 
// — > X can be adapted to other symmetries.) Appropriate random-matrix ensembles 
follow by choosing F = F^ from the standard circular orthogonal or unitary ensembles 
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(COE or CUE), respectively [T7JEB]- We denote these circular ensembles with VT and 
VTT' symmetry as CShSxE and CSrSxE', respectively 



2-4- Overview of characteristic parameters and scales 

In summary, each RMT ensemble is specified by the symmetry S-rSx, as well the 
following 4 dimensionless parameters: the number of modes M — 1/ At in each of the 
two parts of the system, the relative size a = N/M of the interface, the transparency T 
of the interface (encoded in 7 or 7), and the degree of non-hermiticity n/E^, where 
Et = NT/2ttA is the Thouless parameter mentioned in the introduction. In the 
Hamiltonian variants of RMT, A = tt/M and = Ta/2, while in the quantum 
map version with r = 1, A = 2ir/M and = Tot. 



3. Eigenvalue distribution in the large M limit 

In order to get insight into the distribution of eigenvalues in the complex plane, and 
the conditions under which they may accumulate on the real axis, we first consider the 
limit of a large number of internal modes M — > 00, at fixed a = N/M, h/Et and T. 
For an optical system, this limit is realized by decreasing the wavelength (increasing the 
frequency) in a given resonator geometry while keeping the absorption and amplification 
rate \i at a wavelength-independent value. In random-matrix theory, this limit can be 
approached via systematic diagrammatic expansions, where the leading order captures 
the averaged eigenvalues density neglecting fluctuations on the scale of the level spacing 
[25| [29l I3TH [31] ■ We now adapt this approach to the symmetries in question. 



3.1. Generalized Pastur equation 

The effective Hamiltonian ri generally possesses complex eigenvalues, and the complex- 
analysis nature of the method employed below suggests to denote these as z. The 
distribution of eigenvalues in the complex plane can then be written as 

1 1 dtrg n 

p{Z > Z) = 2Mn^z^' (8) 
where Qu denotes the 2M x 2M-dimensional top-left block of the 4M x 4M-dimensional 
matrix Green function 

-1 

(9) 

The limit A — > is implied to be taken at the end of the calculation. 

In order to work out the random-matrix average we expand the matrix Green 
function as a geometric series 

00 

G = U- 1 ^(--Ho^r 1 )™, n = diag(#, H, H, H), (10) 

n=0 
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o 

iA 

-r 
z* + i/i y 

where the momentarily stipulated form of Hq holds for the GOOE and GUOE' ensembles 
(the other ensembles are discussed thereafter). 

The average can now be carried out by contractions of the Gaussian random 
variables in H, which can be represented diagrammatically. The leading order (the 
planar limit) is given by rainbow diagrams in which the contraction lines do not cross, 

g = u- 1 +u- 1 n gn g, (12) 



V 



(14) 



which sum up to 

Q = U- X + U-\G ®\ M \G ^ G = {U- 1 -G®l M )-\ (13) 

where G = jjtr mG is a reduced 4x4 matrix Green function. 

The matrix U~ l — G <g> 1m on the right-hand side of (fi~3l) separates into N blocks 
of the form (w 7 — G) and M — N blocks of the form (uq — G), where 

/ z + i/i —7 iA \ 

—7 z — i/i iA 

iA z* — i/i -7 

iA -7 z* + i/i / 

We thus can invert each block separately, and take the partial trace on both sides. This 
leads to the generalized Pastur equation 

G = a(i h - G)- 1 + (1 - a)(u - G)' 1 (15) 

for the GOOE and GUOE' ensembles. 

For the GUOE ensemble, (ED holds with n = diag(If, H*, H, H*) = 
diag(i/, H T , H, H T ). The transpositions reduce the number of rainbow diagrams in 
the Gaussian average, which leads to the modified equation 

G = a(n 7 - P1GP1 - P2GP2)- 1 + (1 - a)(u - P 1 GP 1 - P 2 GP 2 )-\ (16) 

where P l = diag (1, 0, 1, 0) and P 2 = diag (0, 1, 0, 1). 

For the GOAE ensemble, the random matrix A has to be incorporated into 
%q = diag(if + A, H + A, H — A, H — A) , while U is replaced by Li = W| M= o- Instead, 
/i appears via the contractions of A. This leads to the condition 

G = a(u^ -G + RGR)- 1 + (1 - a)(u - G + RGR)' 1 , 

where R = /idiag(l, 1, —1, —1) and m 7 = ^7^=0- 

Finally, in the GOAE' ensemble we have Hq = diag(if + A, H - A, H - A, H + A), 

and 

G = a(w 7 -G + R'GR'Y 1 + (l — a)(u — G + R'GR')~\ (18) 
where R' = /i diag(l, —1, —1,1). 



(17) 
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3.2. Solution of the generalized Pastur equation 



The condition (fl5|) can be rephrased as 

(w 7 — G)G(uq — G) = it(i_ a ) 7 — G, 



(19) 



thus, a third-degree matrix polynomial, and the versions (TT6l)-(jT8l) can be rewritten 
analogously. If we write G in terms of its 16 components and eliminate these successively, 
we end up with a polynomial of a very large degree, which prohibits an exact analytical 
solution. Therefore, we pursue a semi-analytical approach which starts with the exact 
solution G in the uncoupled case a = 0, where Go(uo — Go) = I4 (independently of 
7). This is solved by G = u /2 + (mq/4 — I4) 1 / 2 , where the square root of the matrix 
K = Mq/4 — 1 4 = Vdiag k n V~ l is defined by diagonalization, K 1 ' 2 = Udiag {^szy/h^V^ 1 , 
which here can be carried out explicitly as K decouples into two independent 2x2- 
dimensional blocks. The correct branch is determined by the limit Gq(z = z* = 0; A = 
0) = —ia x <g> 1 2 . 

To describe the following steps let us denote the desired solution of (fT5T) as 
G(z,z*;a, X). Now, we proceed as follows: (i) We determine G(zq, z$; 0, A) = 
Gq(zo, Zq] A) for a fixed value of zq (e.g., z = 0) and a finite value of A (concretely 
chosen to equal the eventual value of a, as we expect the support of the spectrum to be 
of that order), (ii) For values a increasing in small increments from zero to the desired 
final value, we solve (Tl5|) numerically for G(z, z*; a, X), where the initial condition is 
taken as the solution from the previous step, (iii) Analogously, we next decrease the 
value of A to a small value (here taken as 0.001; keeping A small but finite regularizes 
branch cuts). The same procedure can be applied to solve f|T6|) - f|T8|) . 

In practice, we find that a reliable numerical approximation of the desired solution 
G(z , Zq] a, 0) is obtained in a few (about 10) steps. We can next keep a and A fixed but 
vary z in small steps to sample the complex z plane. Furthermore, by considering z* as 
a formally independent variable we can also obtain the numerical derivatives required 
for the calculation of the eigenvalue probability density 



The lower panels in figure [2] illustrate the resulting eigenvalue distribution for 
the case of the GOOE and GUOE' ensembles, governed by ( fT5l) . as density plots for 
a = 0.2 and T = 1 (7 = 1), for two different values of /1. The top panels show the 
eigenvalues of 20 random GOOE matrices with M = 400, iV = 80. Figure [3] shows 
the corresponding results for the GUOE ensemble, as well as results for the GOAE and 
GOAE' ensembles for a single value of /1, obtained from (fl6|) . (fl7|) and (jTBl . respectively. 
In all cases there is excellent agreement, which includes details such as the branches of 
the eigenvalue support extending along the real axis in the range |Rei£| > 2, observed 
for the GOOE and GOAE' ensembles. (We also confirmed quantitative agreement by 
comparing histograms at fixed KeE.) 

The one feature which is not captured by the diagrammatic expansion is the 
visible accumulation of eigenvalues in the whole range jRe^l < 2 along the real axis. 



p(z, z*) 



J__d_ 

2^dz~ 



-(Gu + £22)- 



(20) 



Random-matrix theory of resonators with VT or VTT' symmetry 



9 




Figure 2. Eigenvalue distributions in the complex plane for the GOOE 
ensemble (representing an optical system without magneto-optical effects and uniform 
amplification or absorption), for N/M = a = 0.2, T = 1 and (a) /i = 0.1 = Et as well 
as (b) fj, = 0.2 = 2£'t- The scatter plots in the top panels are obtained by numerical 
diagonalization of 20 matrices T-L taken from the GOOE with M = 400 and N = 80. 
The results in the lower panels are obtained from the generalized Pastur equation (|15l) , 
as described in the text. These results also apply to the GUOE'. 




Figure 3. (a,b) Same as figure[21 but for the GUOE ensemble. (c,d) Analogous results 
for the GOAE and GOAE' ensembles, respectively, with fi = 2Et. 
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However, we find numerically that with the present scaling of parameters (/i/£r fixed 
independently of M), the fraction of these eigenvalues amongst all eigenvalues steadily 
decreases oc M _1//2 as M increases, indicating that the real component of the spectrum 
indeed becomes negligible in the limit M — > oo. This is consistent with the earlier 
prediction for the GOOE and the GOUE [19] that the transition to a complex spectrum 
happens for a characteristic value /ipx which is much less than if M is large. These 
features render the transition out of the reach of the described diagrammatic approach. 
In the following sections, we will first study the transition in detail based on direct 
numerical sampling and diagonalization of the random-matrix ensembles, and then 
extend the perturbative treatment of [19] to quantify the dependence of /apt on T, 
M and N. 

4. Transition from a real to a complex spectrum 

The transition from real to complex-conjugate pairs of eigenvalues can be quantified by 
considering the fraction f c of eigenvalues which are complex; f c = indicates a fully 
real spectrum, while / c = 1 if the spectrum is fully complex. We determine this fraction 
numerically as a function of the non-hermiticity parameter /i, while keeping T, M and 
N fixed. In the Gaussian ensembles, the energy levels are taken only from the central 
region KeE pa of the spectrum, where A ~ ir/M is approximately constant. This 
eliminates the spectral edge effects observed in some of the eigenvalue distributions in 
the previous sections, and represents the typical physical conditions met in the short- 
wavelength limit of realistic resonators (where the effective level index is large, and the 
spectrum is not bounded from above). Furthermore, for Sx = O (uniform absorption 
and amplification), we compare the results to the circular ensembles, as in these the 
mean level spacing is energy independent. (For Sx = A the quantum maps are less 
convenient.) 

Our eventual goal is to characterize the transition in the different ensembles by the 
coupling dependence of the critical scale 



which we identify via / c (a*pt) ~ 1/2 (without requiring exact equality). The scale /io is 
chosen such that the function g(T) does not depend on M and iV if M ^> N ^> 1 (with 
possible exceptions for weak coupling T < To,a, as specified below). By varying M and 
iV independently we find that this scale depends on the symmetry Sx = O or A, and 
can be suitably written as 



This being fixed, we set M = 200 and N = 40 and determine f c (n) for various values 
of T with \i measured in units of the appropriate /io, as shown in figures HHTJ and focus 
the discussion on the ensemble-specific form of the scaling function g(T) in ( 12TT) . 



/i PT = g(T)na, 



(21) 




VNA/2n = /i , (OO, UO and UO'); 
\/MA/27r = fi A , (OA, OA'). 



(22) 
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In the GOOE and COOE (figure 0]), we find that the transition becomes coupling- 
independent as soon as T > T G , where 



In this regime g(T) ~ 1; the slight T-dependence still observed in the plots are finite-size 
effects, which disappear if M and N are further increased. (However, in this limit Tq 
becomes very small, so that the behaviour for T < To would be difficult to illustrate; 
the chosen values of M and N thus constitute a suitable compromise.) In the GUOE 
and CUOE (figure EJ), on the other hand, the transition displays coupling dependence 
throughout the whole range of T. 

These results for 00 and UO symmetry agree with the predictions in [19] . which 
we systemize in the following section to develop a microscopic picture that also applies 
to the remaining ensembles considered here. The numerical results for these cases are 
as follows. 

In the GUOE' and CUOE' (figure H}, the transition displays a similar coupling 
independence as in the GOOE and COOE, with /xpx (and thus g(T)) roughly scaled 
down by a factor of about \/2. 

In the GOAE and GOAE' (figure [7]), the scale /xa applies. In both ensembles, there 
is almost no coupling dependence throughout the whole range of T, with exception of 
the weak-coupling regime T <C Ta in the GOAE', which is now delineated by 



Interestingly, in this regime g(T) decreases with increasing T, which amounts to an 
anomalous behaviour — the coupling between the resonators enhances the fragility of the 
real spectrum, in contrast to the usual situation where increasing the coupling furthers 
the balance of the non-hermitian effects in the system. 

Note that in the limit M = N/a — > oo studied in section [3], Ho,a/Et -)■ as well 
as T ,a 0. 

5. Relation to level spacing statistics 

To explain the observations of the previous section, we now develop a microscopic picture 
of the transition from a real to complex spectrum. This is based on two ingredients, the 
level statistics in the hermitian limit fi = 0, and the interaction of the eigenvalues on 
the real axis as /i is increased. Focussing on these two ingredients is motivated by the 
fact that the formation of complex eigenvalues requires two real eigenvalues to coalesce. 
The required degree of nonhermiticity thus depends on the distance of the levels at 
H = 0, and the typical size of the matrix elements which mix the levels. The ensembles 
studied here display different degrees of level repulsion and level mixing, which also 
depend on the coupling strength T, and our aim is to show that these characteristics 
are consistent with the coupling dependence of the complex fraction f c reported in the 
previous section. 



T ~ l/N. 



(23) 



T A ~ l/N 2 . 



(24) 
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(a) GOOE (b) COOE 




Figure 4. Average fraction f c of complex eigenvalues as a function of the absorption 
rate fx, scaled to /io = V r NA/2n. Panel (a) shows results for the GOOE around 
KeE = 0, while panel (b) shows results for the COOE, obtained in both cases by 
numerical sampling of the ensembles with M = 200 and N = 40. The different curves 
correspond to different transparencies T of the interface between the amplifying and 
absorbing resonators. For T > Tq ~ l/N, the initial stages of the transition to a 
complex spectrum is coupling-independent. 



(a) GUOE (b) CUOE 




Figure 5. Same as figure [H but for the GUOE (a) and CUOE (b). Here a clear 
dependence on T persists throughout the entire range of this parameter. 

We start with some preliminary observations that justify to separate the problem 
of level spacing statistics at jj, — from the problem of level mixing at finite \i. First, 
we note that at fi = 0, T = 0, the system consists of two uncoupled passive resonators, 
which both have an identical real spectrum. In order to inspect how this degeneracy 
is lifted, we pass over to a "P-symmetric basis, "Hp = UHU, where U = 2~ l l 2 (a x + a z ) 
diagonalizes cr x , which gives 

Uv=\ H + T "^t, I (GOOE and GUOE') (25) 

\ —ifi h — r i 
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(b) CUOE' 



M/M-o 




Figure 6. Same as figure 0] but for the GUOE' (left panel) and CUOE' (right panel). 
Again, a coupling-independent regime emerges for T > Tq ~ 1/N. 




(b) GOAE' 
1 



0.5 
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Figure 7. Same as figure 0J but for the GOAE (a) and the GOAE' (b), and with fx 
now given in units of lla — \/~MA/2ir. The transition is almost coupling independent, 
with exception of the weak-coupling regime T < Ta ~ 1/N 2 of the GOAE', where the 
characteristic scale /zpt increases as T — >• 0. 



Hr I H + ^ + T „ , r ) (GOAE) (27) 






^= ^" A r F I r I (GOAE') (28) 

In the hermitian limit ji = (implying also A = 0), all these transformed 
Hamiltonians are block diagonal, with exception of (126]) for the GUOE. This is the 
case because in the other cases V is an exact symmetry; moreover, for the GOOE, 
GOAE, and GOAE', T and T' hold separately if \i = or A = 0. These properties lead 
to different level statistics in the hermitian limit, which in all cases but for the GUOE 
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Figure 8. Level spacing distribution of real eigenvalues in the hermitian limit of 
the GOOE (curves) and COOE (points), obtained by numerical sampling of these 
ensembles with M = 200 and N = 40. These results also apply to the GOAE' and 
COAE'. 



involve the superposition of two non- interacting level sequences E + and E~ , obtained 
from H + T and H — T, respectively. The two sequences are degenerate at T = 0, but 
are modified by T, which perturbs the two sequences in different ways, and because of 
its positive definiteness also induces an approximately rigid shift. 

In the non-hermitian case (finite /i or A) the V symmetry remains exact in the 
GOAE, while T' symmetry remains exact in the GOOE and T symmetry remains exact 
in the GOAE and in the GOAE'. These differences are reflected in the matrix elements 
that mix the level sequences. For this, we recall that in almost-degenerate perturbation 
theory, the effective Hamiltonian of two adjacent levels Ei and Ej is 



where V is a generic perturbation. The perturbation theory is straightforward at small 
T, but as T increases levels display exact or avoided crossings. One can then still base 
estimates by stipulating a typical spacing A of two adjacent levels at finite T, and 
restricting the perturbative analysis to the wavefunction overlap |19j . 

As a backdrop for the ensemble-dependent discussion of the details, we show in 
figures [8] and [9] numerically evaluated level-spacing statistics -P(s) at ji = 0, where s is 
the distance between adjacent levels. In the Gaussian ensembles, we focus again on the 
bulk of the spectrum (close ReE rs 0); for Sx = O these results are also compared with 
results from the circular ensembles (as before M = 200 and N = 40). 

In figure El P(s) is shown for the GOOE and COOE. The data applies to the full 
spectrum of "H, thus, the superposition of the sequences E + and E~ of H-p at /x = 
(with mean level spacing s = A/2), which is appropriate as these sequences become 
mixed by finite /i. At small coupling T < Tq ~ 1/N, the statistics is dominated 
by the closeness of levels which degenerate at T = 0. Based on (1291) . one then finds 
NVTA/2tt, thus g(T) ~ y/lVT < 1. For T > T , levels in the two sequences 




(29) 
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(a) GUOEandCUOE (b) GUOE ' and CUOE ' (c) GO AE (one block) 




s/A s/A s/A 

Figure 9. Same as figure [SJ but for the GUOE and CUOE (left panel), the GUOE' 
and CUOE' (middle panel), and a single block H + T, corresponding to the levels that 
are mixed by finite A in the GOAE (right panel) 



cross (giving P(0) ~ 1/A), and the spacing statistics quickly converges to a coupling- 
independent form. At finite /i, adjacent levels of the different sequences E£ and 
with \E+ — E m \ ~ A are mixed by a matrix element of squared size \ti{ipn\' l f J m)\ 2 ~ A* 2 /-^" 
|19j . which becomes comparable to A 2 at /ipx ~ \/NA/2ir = /xq, up to factors of order 
unity because of the influence of the level fluctuations. This qualitatively explains the 
approximate coupling-independence of / c (/i), observed for this ensemble in the previous 
section (figure HJ). One can interpolate between the weak-coupling and finite-coupling 
regimes by setting g(T) « (1 + l/NT)- 1 ' 2 . 

The left panel of figure [9] shows the analogous result for the GUOE and CUOE. 
Here, finite V also induces a mixing of the originally degenerate sequences E + and E~ , 
which results in a coupling-dependent level repulsion (with P(0) = 0). This corresponds 
well to the persistent coupling dependence of / c (/i), observed in figure [5j Based on (j29|) . 
one now finds /ipT ~ v / A 7 TA/2tt P2], thus g(T) ~ v 7 ^, which holds across the whole 
range of T, up to modifications of order unity, which we now can associate to the 
influence of the level statistics. 

In the GUOE' and CUOE' (middle panel of figure |H]), on the other hand, levels 
are again not mixed by finite T. Thus, a coupling-independent level statistics again 
emerges for T > To, with is similar to the result for the GOOE/COOE (figure EJ, 
with -P(O) = 1/A. The modal value is shifted to slightly larger s, in accordance with 
the larger degree of level repulsion in the standard GUE/CUE [T7]. This behaviour 
corresponds well to the approximate coupling-independence of / c (/i) in figure |6j We 
find in the perturbative treatment that g(T) is the same as in the GOOE, up to a 
possible factor of order unity due to the small differences in the level spacing statistics. 

In the GOAE, the superimposed level sequences display the same statistics as in 
the GOOE. However, A does not mix these sequences; instead, eigenvalue coalescence 
must happen within a given sequence. Therefore, we consider the spacing within a fixed 
sequence, which is shown in the right panel of figure |9] (here the mean levels spacing is 
s — A). The result is almost indistinguishable from the standard Wigner distribution 
of the GOE (with P(s) oc s for small s) [TTJ [18], as the main effect of T is a rigid shift; 
small deviations appear only for T«l. This agrees well the corresponding behaviour of 
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/ c (aO in figure [7(a). Perturbatively, the mixing is given by overlaps a mn = 

of size a^ nn = n 2 /M, which must be of order A 2 for eigenvalues to coalesce. Thus, we 

can write /ipt ~ y/MA/2n = fix (up to a factor of order unity), which is independent 



The data in figure M also applies to the GOAE', which shows a similar coupling 
independence, figure [TJb), as the GOOE and GUOE'. The anomalous behaviour at 
small T can be understood from the fact that in this regime eigenvalue coalescence 
predominantly occurs between levels E+, E~ of the two sequences that are degenerate 
at T = A = 0. Thus, perturbatively their eigenvectors ip£ ~ 4>n are identical, and the 
first-order coupling (i/;~\A\^) = because A is antisymmetric (and tp^ is real). For 
T > Ta, on the other hand, the coalescence is between originally non-degenerate levels 
of the two sequences, and /ipx ~ \/MA/2n = fi^ (up to a factor of order unity, and 
again independent of N), as in the GOAE. 

6. Conclusions 

In summary, motivated by recent optical realizations of non-hermitian PT-symmetric 
quantum systems, we identified a number of symmetry classes which can be realized 
in optical resonators and differ by a choice between two generalized time reversal 
symmetries (VT or VTT'), as well as the absence or presence of magneto-optical effects 
in the hermitian and nonhermitian parts of the dynamics. Specifically we considered 
five scenarios, with symmetries termed 00, UO, UO', OA and OA'. 

Our analytical results reveal that in the short-wave limit, the fraction of real 
eigenvalues among all eigenvalues in the spectrum decays to zero at any classically finite 
amplification and absorption rate \i. Based on numerical results, and an extension 
of the perturbative results in [19J, we find that the amplification and absorption rate 
/xpt ~ Hog(T) at which real eigenvalues turn complex is indeed characterized by a scale 



which is classically small but microscopically large. Furthermore, the scenarios differ in 
the dependence g(T) on the coupling strength T between the absorbing and amplifying 
components, which can be explained in terms of the level spacing statistics in the 
hermitian limit fi = 0, and distinct mechanisms of how these levels are then mixed 
when /i is finite. For UO symmetry, the transition is T-dependent over the whole range 
of this parameter. In the 00 and UO' classes, a significant dependence is only observed 
for T < To ~ l/N, while the OA' symmetry class displays an anomalous dependence 
of the transition on the coupling strength in the range T < Ta ~ 1/N 2 . In the OA 
symmetry class, there is negligible coupling dependence over the whole range of T. 

The introduced models possess a structure that respects the constraints imposed 
by characteristic energy and time scales, the physical nature of the amplification and 
absorption, and the accessible coupling strengths of a realistic (possibly semitransparent) 



of N. 




VNA/2tt = fio, (00, UO and UO'); 
y/MA/2n = /i A , (OA, OA'), 



(30) 
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interface. The classification of these models can be straightforwardly extended to 
include any symplectic, chiral, particle-hole like, or additional geometric symmetries, 
as previously discussed in hermitian situations [251 EEl E7J [32j |33l EI] . 

In this work we focussed on ensemble-specific spectral properties. However, as is 
generally the case in random-matrix theory, there are many quantities that are less 
ensemble-specific and should display a large degree of universality. The prime example 
is the spectral statistics of the complex eigenvalues in the bulk of the spectrum. If one is 
interested in such statistics, simpler models can be employed. For example, in scattering 
theory [31] the effective Hamiltonian has a semidefinite antihermitian part, but the bulk 
spectral statistics can be studied via the Ginibre ensemble with complex entries [35J. 
Analogously, the spectral constraints of VT or VTT' symmetry are also obeyed by the 
real Ginibre ensemble, which has a much simpler matrix structure than our ensembles, 
and for which detailed rigorous results are available [361 EH EE] . That this ensemble is 
a good model for other non-hermitian ensembles has been demonstrated, e.g., for the 
case of lattice QCD in |39j. Notably, in this ensemble, in the stipulated limit with fixed 
classical scales, the fraction of real eigenvalues decays as 1/a/M [40J. 

An important constraint in the applicability range of any random-matrix ensemble 
is the requirement that many modes are well mixed by multiple scattering. This is not 
the case in (quasi) one dimensional "PT-symmetric disordered systems, where states 
are localized and the transition happens at much smaller values of [i [4"T1 l4*2l |4"3] . 
Furthermore, even in wave-chaotic settings the multiple scattering can be suppressed by 
dynamical effects, which can lead to systematic corrections for the density of eigenvalues, 
as observed in [20] for the strongly amplified states in a quantum-chaotic system. 
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